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Abstract 

We argue that the random-matrix like energy spectra found in pseudointe- 

grable billiards with pointlike scatterers are related to the quantum violation 

of scale invariance of classical analogue system. It is shown that the behavior 

of the running coupling constant explains the key characteristics of the level 

statistics of pseudointegrable billiards. 
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The concepts of scale anomaly and the asymptotic freedom are among the key features 
of the gauge field theories which describe the interaction of the elementary particles. It is 
less widely recognized, however, that the scale anomaly can be found in a vastly simpler 
setting of one particle quantum mechanics. A particle scattered off a pointlike scatterer in 
two spacial dimension is known to have energy dependent s-wave phase shift defying the 
scale invariance of its classical analogue 0. There exists a sister problem of particle mo- 
tion confined in a hard-wall boundary with a pointlike scatterer inside. When the shape 
of the boundary is a rectangle, the problem belongs to a larger category of systems known 
as pseudointegrable billiards [@-|5|]. This system is known for puzzling statistical properties 
of its energy eigenvalues @,[|]. It is shown through numerical experiments that the level 
statistics of the pseudointegrable billiard resembles to that of random-matrix ensembles || 
which is generally associated with chaotic dynamics ]7). This is in seeming contradiction 
with the absence of chaotic dynamics in classical analogue system. Further, when the levels 
are collected at higher energy region, the level statistics moves toward the Poisson distribu- 
tion which characterizes the integrable classical dynamics. Also, the system tends to show 
more of Wigner-like statistics when the genus of the billiard is increased, that is, in the 
present context, when the number of the singular scatterers is increased. These facts have 
never received sufficient explanations, in spite of several attempted studies based on the 
semiclassical periodic orbit quantization theory P|§. 

In this paper, we argue that the behavior of spectra of the pseudointegrable billiard with 
pointlike scatterers is a direct result of the scale anomaly of the system. Specifically, the 
dependence of the level statistics on the energy and the number of the scatterers is shown 
to be well explained by the high energy behavior of the effective coupling strength of the 
pointlike scatterer. 

We consider a quantum particle of unit mass moving freely inside a boundary B in 
two spacial dimension on which its wave functions are assumed to vanish. We denote the 
eigenvalue and eigenfunction of the system as e n and (f> n , namely 
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- ^V 2 (j) n (x) = e n <j) n (x) (la) 

with 

4>n{xB) = where x B G B. (lb) 
Assuming (j> n {x) to be normalized to unity, the Green's function is given by 

G(°\x,x';u J ) = jr M * )MSf) . (2) 

n=l U ~ £ n 

When the shape of the boundary B is such that the classical motion of the particle is regular 
(as in the case of B being rectangular), the nearest neighbor spacing s n = e n+ i — e n is 
distributed according to the Poisson distribution P(s) = exp(— s). We now place a pointlike 
scatterer at xq. Naively, one defines the scatterer in terms of the Dirac's delta function in 
two dimension 

V{x) = v 5(x — x ). (3) 

Under the scale transformation x — > ax, the potential is transformed as V(x) — > V(x)/a 2 . 
Since this behavior is identical to that of the Laplacian operator in eq. (|TaD , the system is 
scale invariant. One expects, therefore, that the dynamical properties of the system should 
not depend on the energy. Formally, the transition matrix T (T-matrix) in the presence of 
the scatterer V is given by 

T — V + V G (0) T. (4) 

The poles of T give the eigenvalues of the system. Because of the separability of the delta 
potential, (</>„ V 4> m ) = v 0„(xo) <fim(xo), the T-matrix is also separable: 

(<p n T <f> m ) = t(u) (f) n (x ) (pm(x ). (5) 

Apart from the trivial solution u = e n for the case of <f) n (xo) = 0, the poles of T are formally 
given by the roots of the equation 
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with 



_i_ = i_ G (°)(^) = (6) 
t{uj) v y 1 w 



G(°)M = G( )(f ,f ^) = f:^^- (7) 

n=l UJ ~ £ n 



However, eq. (H), as it stands, is meaningless since 



„=i u-e n n=1 uj-e n 

roc 1 
\2\ ' • L 



« (0(x o ) 2 ) po / efe ► oo (8) 

Jo lu — e 

where (4>(xq) 2 ) is the average value of <j) n (xo) 2 among various n. The divergence is brought 
about because the density of states is constant (which we denote po) with respect to the 
energy. To handle the divergence, a scheme for regularization and renormalization is called 
for. The most mathematically satisfying scheme is given by the self-adjoint extension theory 
of functional analysis [|K],|ll[]. Here we just quote the result. After the self-adjoint extension, 
the transition matrix t(u) is given by 

1 (uj - iA) 



t(cu) 1 - e 



dx G (0) (x, x ; u)G (0) (x, x ; iA) 



+Y^^ J dxG^(x,x ;u)G^(x,x Q ;-iA). (9) 

Here, A is an arbitrary scale of the regularization, and (0 < < 2tt) is the parameter 
of self-adjoint extension. With the straightforward calculation, we find that the energy 
eigenvalues of the system - the poles of t{uS) - are determined by the equation 

- - G(u) = (10) 

where 
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n=l 



+ 



U - E n el+A 2 



(11) 



is the regularized version of G(°\uj), and 
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A sin 6 



E 



4>n {X f 



l-oose^ +A 2 



(12) 



is the formal (or renormalized) coupling strength of the scatterer. We stress that in spite 
of purely mathematical construction of eqs. ([10| ) - ([T2|) , it does correspond to the physical 
small-size limit of the problem of a finite-size obstacle. 



Since the series of eq. (11) is convergent, the problem is now well defined. We look at 
the behavior of eq. (|10|) at high energy region uj » A. For a given value of u, we can 
approximate eq. (|TT|) by truncating the summation at n — n x (u>) 

1 



g( w ; 

with an error given by 
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n x (ui) 

E M^o) 2 

n=l 



uj — e 7 



+ 



(13) 



E M^o) 

n=n x 

poo 

(0(x o ) 2 )po / de 

Je x (ui) 

n x (uj) 



1 

+ 



uj — e r 
1 



el + A 2 



+ 



uj - e e 2 + A 2 



(14) 



where we have used £ x {uj) = n x (uj)/p . Therefore, we can set 



n x {uj) = oluj 



(15) 



where a is a constant inversely proportional to the allowable error SG. Once the summation 
is truncated at finite terms, we can rewrite eq. fllTf) as 
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v eff (uj,v) 



E 

71=1 



UJ — E r 







(16) 



with the effective strength v e ff(u,v) defined through 
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n x (uj) 



V eff {UJ,V) 



E M%o) 



n=l 



+ A 2 ' 



(17) 



Comparing eq.(Ti6|) and eq. (^|), one realizes that the problem is now turned into an eigenvalue 
problem with finite basis states <f> n , n = 1, • • • ,n x (uj). Although the system originally has 
no inherent scale, the effective coupling strength v e ff has an energy scale A, and, as a result, 



it acquires energy dependence. This is possibly the simplest example of the scale anomaly 
|Ij]. The effective coupling v e ff is also referred to as running coupling because of its energy 
dependence. When v e g is large for the energy region of the interest, it induces the mixing 
among the basis states, and results in the Wigner form of the nearest neighbor spacing 
distribution P(s) = ^tts exp(— \t^s 2 ). It is known that v e g is large only in one energy region 
determined by the value of v and A [||]. Replacing the summation in eq. ( 17|) with the 
integral as in eq. fll4|) , we have 

v eff (uj,v)tt (18) 

1 - t;(0(fo) 2 )po log y/l + (n x (u)/p A) 2 



At the limit u> — > oo, we have log Jl + (n x (uj) / p A) 2 ps logn x (u>) ~ logcu. We arrive at 

v ef f(u,v) pa — J" — (u; — >■ oo). (19) 

ein ; (0(f o ) 2 )p O logu; 

Remarkably, the formal strength disappears from the expression of the effective strength 
v e ff in the high energy limit. It is now easy to see the reason for the level statistics of 
pseudointegrable billiards becoming more Poisson-like at higher energy region, irrespective 
to the choice of the formal coupling v. The fact that the strength v e ff disappears at the limit 
uj — > oo goes along well with our intuition that, at the classical limit, a pointlike obstacle has 
no effect on the motion of a particle. At this limit, all the wave functions are unperturbed 
and the scale invariance is restored in a trivial manner. We note that when we replace the 
sum with the integral in obtaining eq. (plf), we implicitly assume that the size of the billiard 
is far larger than the scale in discussion. That is the reason why the size of the billiard 
boundary, which obviously breaks the scale invariance, does not appear in our arguments. 

We next consider the case of two pointlike scatterers. We place the scatterers at x and 
X\ with the formal strengths vq and v±, respectively. The eigenvalues of this system are 
determined by 



— -G 00 (u) -G 01 (w) 
-GioM — -Gn (a;) 
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(20) 



where Gij(u) and Gij(u) are defined as 



and 



n=l 



1 e r 

+ 



(21) 



Gy(o;) = g ^(^)0n(^) _ (22) 



n=l ^ ~ £ n 



Let us suppose that two scatterers are placed closely to each other. We truncate the sum of 
Goo(^) an d Gn(u) at n = n x {uj) as before. As for Gq\(uj), which is finite as long as xq ^ x±, 
the same truncation is possible if two scatterers are apart by 



\xi - x \ ~ -— = J P ° , (23) 



since the contributions from higher n, which probe finer length scale than \x\ — xq\, cancel 
among themselves. Bellow n < n x (u), we can approximate 4> n (x ) ~ n (xi), since 4> n (x) 
is slowly oscillating in the distance \x\ — x \. The matrix equation, eq. (PUD, then can be 



reduced to 



with the effective coupling v^(u,v ,Vi) given by 



V e ff[UiV ,V 1 ) n=1 

(2) 
eff 



vyiu^o,^) =v eff (uj,v ) +v eff (u,v 1 ). (25) 

This equation reveals an interesting feature of the system with two pointlike scatterers. If 
Vq differs from v% appreciably, v e ff(u, Vq) and v e ff(u, V\) become large at different energies. 
This means that the particle moving in the billiard cannot see the two scatterers at the same 
time for any given energy. On the other hand, the two scatterers disturb the particle in a 
coherent manner when Vq ~ V\. With eq. (|T9|) , we obtain in the limit of V\ = Vq 



v^(uj,Vo,v ) = 2v eff (u,v ) 



2 

(cu -> oo). (26) 



(4>(x ) 2 }po \ogu 
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Namely, two closely placed pointlike scatterers with the same formal strength act as a single 
scatterer of twice the effective strength. 

If we remove our assumption of two pointlike scatterers located closely, we have to 
deal with eq.(pO|) directly. However, we do not expect an essential change of the character 
of level statistics, since the statistical measures are known to be rather insensitive to the 
precise location of the scatterers. This allows us to generalize our findings to the case of 
more than two scatterers. Let us consider a billiard with an arbitrary (finite) number of 
pointlike scatterers inside. We classify the scatterers according to the magnitude of the 
formal strength; we collect the scatterers with the same order of magnitude of the formal 
strength as a single group. We can expect that the scatterers belonging to one of such groups 
disturb the particle motion in a coherent manner in the energy region determined by their 
own formal strength, while their influence never appears at different energies. In particular, 
in the case of scatterers with a common formal strength, their effects are additive. This 
essentially explains the known behavior of pseudointegrable billiards that by increasing the 
number of singular points one obtains more Wigner-like level statistics. 

We now place our findings in a broader context. The standard approach in quantum 
chaology has been the semiclassical periodic orbit theory P~2| , |T3[ which is a truncation of 
WKB approximation. It cannot be easily applied to the problem with singular scatterers, 
since the classical orbit hitting the pointlike scatterer occupies a set of measure zero in 
the phase space. Our approach, on the contrary, is fully quantum mechanical without any 
resort to semiclassical approximations. Although the present approach is applicable only 
to the system with pointlike scatterers, our results certainly give an insight to the problem 
of quantum level statistics of generic pseudointegrable system. Specifically, we believe that 
it establishes a firm ground to the earlier intuitive notion of "wave chaos" @,f|, that is 
the chaotic motion generated by the uncertainty principle in the motion of wave-natured 
quantum particles. 

The system studied here clearly is a pedagogical example of the quantum violation of 
classical scale invariance. We believe that it enhances the intuitive understanding of the 
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scale invariance and asymptotic freedom which have been viewed as phenomena found only 
in the esoteric theories of elementary particles. Finally, we stress that, even apart from the 
arguments on quantum chaos, we could have arrived at the level statistics in search of the 
measurable consequence of scale anomaly of the system, since it gives a handy way to see 
the effective strength of the pointlike scatterer in the absence of clear-cut measures to gauge 
its effect on the individual states. 

In summary, we have derived the energy and number-of-scatterer dependence of the 
effective coupling strength of pointlike scatterers in billiard systems. We have shown that 
it explains the key characteristics of the level statistics of the quantum pseudointegrable 
billiards. 
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